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NAVIGATION FLIGHT OPERATIONS FOR MARS 

PATHFINDER 

Robin M. Vaughan + , Pieter H. Kallemeyn, Jr. 1 , David A. Spencer* , 

and Robert D. Braun* 

On July 4, 1997, Mars Pathfinder became the first spacecraft to 
land on the surface of Mars in 21 years. Pathfinder was launched 
on December 4, 1996 and spent seven months en route to the red 
planet. This report discusses the navigation flight experience for 
the Mars Pathfinder interplanetary cruise. In particular, orbit 
determination and maneuver design and execution results are 
presented. Special emphasis is given to the navigation role in the 
days and hours leading up to and including the Entry, Descent, and 
Landing (EDL) phase. 

INTRODUCTION AND MISSION OVERVIEW 

Mars Pathfinder was one of the first two missions to be launched under NASA’s new 
Discovery Program. The mission was designed to demonstrate a low-cost, reliable 
system for entering the Martian atmosphere and placing a lander safely on the surface of 
Mars. Pathfinder also carried a science payload that returned data on Mars’ atmosphere, 
meteorology, geology and morphology, and the elemental composition of rocks and soil 
near the lander. A key element of this payload was the Sojourner microrover which was 
the first mobile vehicle ever operated on the surface of another planet. 

The spacecraft traveled along the transfer trajectory from Earth to Mars shown in 
Figure 1. During the seven-month journey between planets, the lander and its science 
payload were enclosed in an atmospheric deceleration module (backshell and heat shield) 
that was attached to a cruise stage as shown in Figure 2. The cruise stage contained the 
propulsion unit, star and sun sensors, solar arrays, and other subsystems needed to 
travel from Earth to Mars. The atmospheric deceleration module consisted of the 
backshell, heat shield, braking rockets, and parachute. The lander contained the flight 
computer, X-band radio system, and all engineering subsystems which were used during 
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surface operations. Telecommunications and navigation tracking during cruise were 
performed using an X-band radio system and medium gain antenna on the cruise stage. 
Navigation tracking data was limited to coherent, 2-way (X-band) Doppler and ranging. 

The Pathfinder spacecraft entered the Martian atmosphere directly from its 
interplanetary trajectory. Initial deceleration in the atmosphere was accomplished by 
aerodynamic braking using the heatshield, followed by deployment of the parachute. 
Small braking rockets were fired at a specified altitude to null the remaining vertical 
velocity. Airbags were inflated surrounding the lander to cushion the surface impact, the 
lander was separated from the backshell and parachute, and allowed to free fall the rest 
of the way to the surface. The lander bounced at least 15 times before coming to rest on 
the surface. The airbags were then deflated and retracted, and the lander was configured 
to expose the payload and rover for surface operations. 

There were 3 requirements on the navigation system for the Pathfinder mission. Two 
of these involved targeting the spacecraft for safe atmospheric entry from the 
interplanetary trajectory and landing at the desired location on the surface. 1,2 The third 
requirement was a limit on the total AV that could be expended for trajectory correction 
maneuvers (TCMs) over the course of the mission. A previous paper 2 described the 
navigation system designed to meet these requirements and presented pre-launch 
analyses of system performance. This paper documents actual navigation performance 
for the mission. The first section is devoted to spacecraft orbit determination (OD) 
topics including force modeling, tracking data quality, and filter configuration. CD 
solutions for the designs of the 4 cruise TCMs are also presented. This is followed by a 
discussion on maneuver design computations and maneuver implementation modes. 
Trajectory targets and the computed velocity change vectors are given for each TCM. The 
design values are compared to actual values from CO solutions reconstructing maneuver 
execution. The third section discusses navigation activities for EDL; results of the final 
CO solutions using tracking data received up to the last minutes before atmospheric 
entry are presented. The final section gives conclusions on the overall Pathfinder 
navigation experience. 


(Note: Actual launch date was December 4, 1996.) 

Figure 1 Mars Pathfinder Interplanetary Figure 2 Mars Pathfinder 
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ORBIT DETERMINATION 


Two important aspects of the orbit determination (OD) task were spacecraft non- 
gravitational force modeling and the quality of the radiometric tracking data. Overall 
force modeling is summarized and details are given on the evolution of the spacecraft 
component model for solar radiation pressure. Variations in Doppler and ranging data 
quality throughout the mission are shown and their impact on the orbit determination 
process is discussed. Modeling of error sources for the Doppler and ranging tracking data 
is described. The filter configuration for in-flight trajectory estimation is compared to 
the pre-launch plans. 1,2 The results of orbit determination solutions for the designs of 
TCMs 1-4 are presented. (OD results from the last 48 hours before entry are covered in 
the section on navigation for EDL.) 

Non-Gravitational Force Models 

Other than the gravitational attraction of the Sun and planets, the only known forces 
acting on the Pathfinder spacecraft during cruise were thrusting events and solar 
radiation pressure. All attitude control turns and spin rate changes were modeled as 
“small forces” - instantaneous velocity changes with magnitudes ranging from 0.5 to 7 
mm/sec. A total of 24 such events were eventually included in the filter runs, covering 
initial turns to reduce the spin rate from 12 to 2 rpm after launch, 17 turns for 
pointing control during cruise, and cruise stage separation just prior to atmospheric 
entry. The four TCMs were modeled as a set of either finite or instantaneous burns, 
depending on their magnitude and execution mode. TCM 1 was modeled as one continuous 
finite burn since it was executed in one segment lasting about 90 minutes. The remaining 
3 TCMs were each executed as multiple burn segments. TCM 2 was modeled as two finite 
burn segments, with the longest segment lasting approximately 5 minutes. TCMs 3 and 4 
were modeled as 3 and 2 instantaneous burn segments, respectively; these were 
relatively small velocity changes with the segments completing in less than 30 seconds. 

Solar radiation pressure contributed a small, but significant, continuous force on the 
spacecraft. Constructing the filter’s solar pressure model proved more difficult than 
anticipated. The shape and reflectance properties of the spacecraft needed to be 
approximated in the orbit determination software using components of various standard 
types such as flat plates and cylinders. A key issue was the modeling for the backshell 
which has a non-standard shape and which was shaded by the cruise stage to a varying 
extent during most of the flight. The challenge was to choose a component type and size 
that captured the contribution of the solar pressure from the illuminated portion of the 
backshell. An additional complication was the inability of the software to model the 
effects of shading due to changing attitude and component configuration. 

Ultimately, a total of four solar pressure component models were used during flight. 3 
All of the models used a flat plate to represent the solar array disk of the cruise stage. 
The launch vehicle adapter and the heat rejection system (HRS) panels, also on the 
cruise stage, were represented by another flat plate and a cylinder in the later models. 
The primary difference in the models was the treatment of the backshell. The first and 
fourth models used a flat plate active for the entire cruise. In contrast, the intermediate 
models invoked the navigation software’s ability to limit the time span over which a 
given component contributes to the solar pressure force. The second model represented 
the backshell as a flat plate active only after March 1997. The third model used two 
separate components for the backshell: a cylinder active before April 1997 and aflat 
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plate active thereafter. These were attempts to account for the changing solar 
illumination on the backshell; the orientation of each component relative to the cruise 
stage components as well as the time of activating components was based on the history of 
the solar aspect angle over the mission. 

Various problems with the filter solutions were encountered when using each of 
these models. Solutions using the early models and estimating reflectance coefficients and 
area scale factors for each component failed to produce a predicted future trajectory 
consistent with the actual tracking data received after the cutoff for those solutions. 
Also, solutions would often produce corrections for the reflectance coefficients that 
would make the new values physically impossible. Experimentation with the filter 
eventually yielded reasonable results with good fits through new data when only the area 
scale factors for each component were estimated. These were the only solar pressure 
parameters included in the OD solutions supporting TCM 4 design and EDL. A more 
serious problem involving switching component types for the backshell was uncovered 
during the last week before entry. Deleting a single pass of Doppler data occurring near 
the April 16th transition date for the backshell component resulted in a significant shift 
in predicted Mars B-plane state. The elimination of these few dozen points over the data 
arc running from February to June should not have made such a large change in the 
solutions. Experiments showed that identical results were obtained if only a single set of 
components were active for the entire data arc, regardless of whether or not that 
Doppler pass was included. Based on this observation, the navigation team opted to model 
the backshell as a flat plate over the entire data arc for all OD solutions following TCM 4 
execution and leading up to EDL. 

Small modeling errors and potential unknown forces (such as outgassing) were 
accommodated in the filter by including a set of stochastic non-gravitational 
accelerations. These parameters were tightly constrained by the choice of a priori 
sigmas so that they were not confused with known forces already in the filter model. 

Radiometric Data Quality 

The quality of the radiometric tracking data was crucial to obtaining the required 
navigation accuracy for Mars Pathfinder. The enhanced filter configuration depended on 
low-noise, high-precision ranging measurements to safely meet the atmospheric entry 
trajectory requirements. Unfortunately, ranging performance was beset with a series of 
problems during the first 6 months of cruise, including unexpected biases and high noise 
levels. In contrast, the Doppler system performed well throughout the flight. The only 
significant influence on Doppler data quality was the signature imposed by spinning the 
spacecraft. 

The Mars Pathfinder spacecraft was spin-stabilized, constantly rotating about an 
axis pointing upward through the center of the aeroshell, roughly normal to the circular 
solar array panel. The polarization of the radio signals and a slight offset of the antenna 
boresight from the spin axis resulted in both a fixed bias and an amplitude modulation 
appearing in the 2-way Doppler data. 4 The Delta launch vehicle released the spacecraft 
with the spin axis pointed towards the sun. The original plan was to orient the spacecraft 
such that the spin axis was very nearly Earth-pointed during the first week after launch 
and to maintain that attitude throughout cruise. The Doppler amplitude modulation is 
negligible for small off-Earth angles and the navigation team intended to ignore it when 
processing the data. However, problems with the sun sensors resulted in the spacecraft 
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being kept in its Sun-pointed attitude until just after the execution of TCM 1 on January 
10, 1997. Thus, for the first 30 days after launch, Doppler amplitude modulation due to 
changing Earth aspect angle was clearly visible in the Doppler data and was used to 
verify the post-launch spin rates and spacecraft attitudes. 

Figure 3 shows the post-fit Doppler residuals over the entire cruise period from 
launch to entry (December 1996 - July 1997). The spin signature due to off-Earth 
pointing is clearly seen in the passes before TCM 1 which have larger scatter, ranging 
from 0.01 to 0.03 Hz (0.18 to 0.53 mm/s). Early OD solutions treated this signature as 
simply an increase in data noise and the data were consequently deweighted. Typical 
Doppler weights for the period before TCM 1 were between 0.10 and 0.025 Hz (1.8 and 
0.44 mm/s), assuming a 60-sec count time. The Doppler passes in Figure 2.2 after 
TCM 1 appear much more compact, signaling the return to an Earth-pointed attitude. The 
scatter of points within each pass decreased to 0.00173 Hz (0.031 mm/s) between 
TCMs 1 and 2. The Doppler noise level continued to decrease, averaging 0.0010 Hz ( 
0.018 mm/s) per pass between TCMs 2 and 3 and 0.0014 Hz (0.025 mm/s) per pass 
after TCM 3. These later Doppler data were ultimately weighted at 0.006 Hz (0.11 
mm/s), again assuming a 60-sec count time. 

The constant bias imparted by the spinning spacecraft has been removed from the 
data before filtering so that it does not appear in Figure 3. At the nominal spin rate of 2 
rpm, the expected value of the bias is 0.0725 Hz. However, the spacecraft attitude 
control system did not maintain the spin rate exactly at this value. A stochastic Doppler 
bias parameter was routinely added to the OD solutions to account for slight differences 
between the actual spin rate and the value assumed for the pre-applied constant bias. The 
time history of this parameter in later OD solutions showed a consistently non-zero and 
predictable trend. This prompted an update to the pre-applied Doppler bias model based 
on the actual history of spin rate values from spacecraft telemetry. The final, critical 
OD solutions for EDL used a linearly varying bias between TCM 2 and 3 and a constant 
bias for a spin rate slightly below 2 rpm thereafter. Constant biases were also applied 
for known spin rate values in the period before TCM 2. 
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Figure 3 Post-fit Doppler residuals from launch to Mars arrival 

Pre-launch requirements levied on the ranging system imposed a la noise (or 
jitter) level no greater than 1 .0 m (~7 Range Units or RU) within each pass. However, 
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the range noise persisted at higher levels until late May as seen in Figure 4 which plots 
the the post-fit ranging residuals for the entire cruise.* The first data returned after 
launch exhibited la noise between 13 and 27 RU (1.86 - 3.8 m). The noise decreased 
after TCM 1 when the spacecraft was oriented to an Earth-pointing attitude. The DSN 
ranging system uses the Doppler signal to assist in forming the ranging measurements 
and the spin signature imposed on the Doppler evidently degraded the early ranging 
measurements. Although the noise level had improved after TCM 1, it remained slightly 
above the desired value, varying between 7.4 and 11.7 RU (1.06 - 1.67 m), until 
March. Following nominal operational plans, the spacecraft’s ranging system 
configuration was changed at this time, significantly lowering the la range noise level 
for passes from DSS 45 (Canberra, Australia) and DSS 65 (Madrid, Spain) to 2.23 RU 
(0.32 m). But the passes from DSS 15 (Goldstone, CA) remained at a high noise level, 
eventually increasing to about 23.3 RU (3.3 m) by the end of May. A concentrated effort 
by DSN personnel eventually traced this problem to a faulty hardware component at DSS 
15. This component was replaced by May 30, and subsequent ranging passes from DSS 
15 had noise levels comparable to the other two DSN sites. The last 35 days of ranging 
data before Mars arrival had an average la per-pass noise of 4.0 RU (0.57 m), well 
within the required level. 

The range data weight in OD solutions was varied throughout the data arc to account 
for the changing noise level. Data before March 7 were typically weighed at 21 RU ( 3 
m). DSS 15 range data continued to be weighed at this larger value up to May 30, while 
DSS 45 and 65 data were weighed more tightly at just 4 RU (0.57 m). After May 30, 
data from all sites was weighed at 4 RU. 
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Figure 4 Post-fit range residuals from launch to Mars arrival 

Ranging biases of various types were routinely employed as part of the OD filtering 
strategy. These included a bias of constant value for all ranging data and stochastic biases 
whose values varied for each pass. Pre-launch expected values for the pass-dependent 


* The residuals in Figure 4 are zero-mean since biases for each pass were estimated in the filter solution. 
Plots of residuals without these biases are presented in Figure 5 as part of the discussion on ranging 
biases. 
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biases were on the order of 10 RU. In-flight OD solutions consistently yielded values 
from 20-50 RU for these biases. The apriori sigma for these parameters in the CD 
filter was increased to 40 RU to reflect these larger values 

Fixed biases for all passes from a specific DSS site were also available in the CD 
software, but it was expected that these site-dependent biases would not be needed. 
However, in the weeks between TCM 1 and 2 it became apparent that the ranging data 
from DSS 15 was offset off from the DSS 45 and 65 data by a rather large amount. This 
is illustrated in Figure 5 which shows post-fit ranging residuals over the mission 
without any estimated biases. A similar situation occurred between TCM 2 and 3 where 
the DSS 45 data appeared biased from the DSS 15 and 65 data (not shown in Figure 5) . 
The causes of these biases were identified and appropriate biase values were applied to 
the affected passes in subsequent CD solutions. A final site-dependent bias was observed 
in the last 35 days of tracking data between the DSS 65 ranging passes and the 15 and 
45 passes. The cause of this last bias was unresolved before landing and a set of site- 
dependent range biases was added to the final CD solutions to account for it. The apriori 
sigma on the DSS 15 and 45 bias parameters was set at 7 RU; while the DSS 6 5 
parameter’s sigma was set at 30 RU. The estimated value of the DSS 65 bias in these 
runs ranged from -15 to -20 RU relative to the DSS 15 and 45 passes. 
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Figure 5 Range residuals from launch to Mars arrival without estimated 

range biases 

Orbit Determination Filter Configuration 

Traditional orbit determination filters have been configured with a mix of (constant) 
bias, stochastic and considered parameters. The set of bias parameters typically includes 
the known non-gravitational forces acting on the spacecraft. Pathfinder’s CD runs 
followed this pattern, estimating solar pressure, velocity changes from TCMs and 
attitude control events, and the spacecraft epoch state as bias parameters. DSN station 
locations were also estimated as bias parameters as were Earth and Mars ephemeris 
parameters in runs performed within 10 days of arrival. The Pathfinder filter 
configuration differed from the traditional approach in that there were no consider 
parameters. Instead, new stochastic parameters were introduced to model data error 
sources that were formerly treated as consider parameters. An important group of these 
new stochastic parameters were the pass-dependent biases for the range and Doppler 
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data. The introduction of these biases for the range data is called the precision range 
technique and allowed the data to be weighted at its inherent accuracy of a few meters. 
The pass-dependent Doppler biases were included to account for small deviations from 
the nominal spin rate. Constant biases were sometimes used in combination with these 
stochastic biases. A constant range bias for spacecraft hardware delays and constant 
biases for each of the DSN sites were used in some runs. Atmospheric media 
transmission effects on the radio signal and Earth timing and polar motion parameters 
were also treated as stochastic parameters. The remaining stochastic parameters were 
components of accelerations acting in the radial and transverse directions. One or all of 
these components were included to account for small, unknown errors in the force model. 
The success of Pathfinder navigation has provided valuable in-flight validation of this 
enhanced filter configuration, with its associated precision ranging technique. 1 ' 2,5,6 

The filter configuration for the OD solutions selected for the final design each of the 4 
TCMs is summarized in Table 1. In general, this proved remarkably similar to the p re- 
launch planned configuration. 2 The two major changes over the course of the mission 
were the set of solar pressure parameters and data biases. OD runs for TCMs 1 and 2 
used earlier solar pressure models having fewer components; area scale factors and 
reflectance coefficients for each component were estimated. Runs for TCMs 3 and 4 used 
later models and estimated only the area scale factors for their components. Pass- 
dependent Doppler biases were used up to the run for TCM 4 design. The new Doppler 
spin signature model was introduced at this time, replacing these stochastic parameters 
with more accurate pre-applied adjustments and a single stochastic Doppler bias 
parameter. The pass-dependent range biases were augmented with site-dependent biases 
to accommodate the observed range data behavior. A similar configuration was adopted 
for subsequent OD runs supporting critical EDL operations; this configuration is 
highlighted by the shaded lines in Table 1. (Representative data weights were presented 
in the preceding sections). 

Orbit Determination Results for TCM Designs 

The results of the OD solutions for design of TCMs 1-4 are summarized in Table 2. 
Predicted spacecraft state at Mars arrival is given as coordinates in a Mars-centered B - 
plane coordinate system. The B-plane is perpendicular to the incoming asymptote of the 
spacecraft’s trajectory (called the S direction) and contains the two axes R and T. T is 
defined as the intersection of the B-plane with the Mars mean equator of date. Predicted 
position in the B-plane is given in Table 2 as a pair of R and T coordinates for the 
location where the asymptote intersects the B-plane. The vector pointing from the origin 
of the R and T axes (Mars center) to the predicted B-plane position is call the miss 
vector or B-vector. Uncertainty in B-plane position is represented by the size and 
orientation of a la error ellipse around the predicted position. The size of these ellipses 
are given as semi-major and semi-minor axis lengths (SMAA & SMIA), while the 
ellipse orientation is given by the angle to the semi-major axis direction measured 
clockwise from the T axis (0). Predicted position in the S direction is given in Table 2 
by the time of closest approach (TCA) and its associated la uncertainty. 
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Table 1 OD FILTER CONFIGURATION FOR TCM DESIGN AND EDL RUNS 



a priori 

Process Noise 

Correlation time 


Filter Parameter 

uncertainty 

(N/A for bias 

(N/A for bias 

Comments 



(la) 

parameters) 

parameters) 


Bias (Constant) 





S/C Epoch State 

position - 1000 km 
velocity - 1.0 km/sec 

N/A 

N/A 

All runs 

Solar Pressure 
Area Scale Factors 

1.5-50% of nominal 
(0.015, 0.05, 0.50) 

N/A 

N/A 

All runs 

Specular Reflectivity 

0.01 16-0.0165 

N/A 

N/A 

TCM 1 & 2 
runs 

Diffuse Reflectivity 

0.02 - 0.10 
(0.0697, 0.0999, 
0.02, 0.1) 

N/A 

N/A 

TCM 1, 2, & 
3 runs 

Earth/Mars ephemeris 

Mars: 1.0 (radial) 

7.0 (transverse) 
7.0 (normal) 

N/A 

N/A 

TCM 4, 
EDL runs 

TCMs 

TCM 1 - 1 finite burn 

pointing - 0.5 deg 
magnitude - 5% 

N/A 

N/A 

TCM 2 runs 

TCM 2 - 2 finite burns 


N/A 

N/A 


TCM 3 - 3 inst burns 

2% in AV components 

N/A 

N/A 

TCM 4 & EDL 
runs 

TCM 4 - 2 inst burns 

3.33% in AV 
components 

N/A 

N/A 

EDL runs 

Attitude Control A Vs 

2.0E-6 km/sec 

N/A 

N/A 

All runs 

S/C Range Bias 

400 - 1000 RU 
( 57 - 143 m ) 

N/A 

N/A 

TCMs 1, 2, 
& 4 and EDL 
runs 

Site Range Biases 

7 RU (1 m) DSS 15,45 
30 RU (4.3 m) DSS 65 

N/A 

N/A 

TCM 4 & EDL 
runs 

Station locations 

10 cm 

N/A 

N/A 

All runs 

Stochastic 





Radial Accels 

1 - 5E-6 mm/sec 2 

1-5E-12 

km/sec 2 

0 (white noise) 

All Runs 

Transverse Accels 

1 - 5E-6 mm/sec 2 

1-5E-12 

km/sec 2 

0 (white noise) 

TCM 2 & 3 
runs 

Pass Range Biases 

40 RU (5.7 m) 

40 RU (5.7 m) 

Variable 

All runs 

Pass Doppler Biases 

0.05 Hz 
(0.89 mm/sec) 

0.05 Hz 
(0.89 mm/sec) 

Variable 

TCMs 1, 2, 3 
runs 

Doppler Bias 
(spin signature) 

0.0007 Hz 
(0.01246 mm/sec) 

0.0007 Hz 

15 days 

TCM 4 & EDL 
runs 

Troposphere 

5 cm 

5 cm 

3 hours 

All runs 

Ionosphere 

4 cm (day) 

1 cm (night) 

4 cm (day) 

1 cm (night) 

5 hours (day) 

24 hours (night) 

All runs 

Polar motion 

10 cm 

10 cm 

48 hours 

All runs 

Earth rotation 

0.001 sec 

0.001 sec 

24 hours 

All runs 





















































































































Table 2 B-PLANE STATE, TIME OF CLOSEST APPROACH, AND THEIR la 
UNCERTAINTIES FOR CRUISE OD SOLUTIONS 


CD 

Data Arc 

B*R 

B*T 

TCA 

Solution 


(km) 

(km) 

(UTC) 

TCM 1 

1 2/9/96-1/2/97 
(24 days) 

279321 .23 

386696.07 

7/3/97 09:14:57 

TCM 2 

1 2/4/96-1/26/97 
(53 days) 

8491 .72 

9767.62 

7/4/97 15:45:26 

TCM 3 

2/4/97-4/29/97 
(84 days) 

-1790.15 

-4494.02 

7/4/97 16:52:52 

TCM 4 

2/4/97-6/23/97 
(139 days) 

-1852.63 

-4250.28 

7/4/97 16:54:21 

CD 

la SMAA 

la SMIA 

e 

la TCA 

Solution 

(km) 

(km) 

(deg) 

(sec) 

TCM 1 

100.64 

11.066 

24.44 

25.93 

TCM 2 

249.03 

234.915 

31 .54 

80.02 

TCM 3 

34.49 

30.717 

106.99 

6.21 

TCM 4 

17.44 

4.192 

109.59 

0.38 


MANEUVER DESIGN AND EXECUTION 


Four trajectory correction maneuvers (TCMs) were planned and executed during 
cruise, the first two within 60 days of launch and the last two within 60 days of Mars 
arrival. This section discusses the choice of target aimpoints, the computation of 
required velocity changes, and implementation modes for each of these maneuvers. 
Maneuver performance is summarized by comparing design velocity changes with CD 
estimates of the actual changes based on tracking data taken during and after the 
maneuvers. The actual total velocity change implemented for the mission is compared to 
the pre-launch budget statistics. 2,7 

Target Aimpoints and Designed Velocity Changes 

The target trajectory conditions for TCMs 1-4 are presented in Table 3 as B*R and 
B»T coordinates in the Mars-centered B-plane coordinate system and an associated time 
of closest approach to Mars. These quantities are chosen for ease of comparison with the 
OD solution results given in Table 2. The launch injection target was biased away from 
the required atmospheric entry aimpoint to minimize the probability of the launch 
vehicle upper stage impacting Mars. The targets for TCMs 1 and 2 were to remove this 
bias, but could themselves be biased from the final aimpoint in order to minimize the 
probability of the spacecraft impacting the Martian surface at a speed greater than 1000 
ft/sec. Analysis of the actual injection accuracy and OD knowledge at the time of TCM 1 
design showed that the impact speed constraint would be satisfied without adding any bias 
to the direct entry target point. This is the reason that the target aimpoint for TCM 1 
listed in Table 3.2 is very close to that of TCMs 3 and 4. Subsequent OD solutions 
reflecting TCM 1 execution showed that this second bias would have to be put back in to 
satisfy the surface impact speed constraint at the time of TCM 2 design, so its target 
parameters differ significantly from those of the other 3 maneuvers. 























































Targeting for TCMs 1, 3 and 4 was based on desired flight-path angle at a specified 
altitude above Mars and landing site latitude and longitude. OD solutions provided the 
estimated arrival geometry and time at the selected Mars atmospheric interface altitude. 
The descent trajectory through the atmosphere then determined the landing site 
associated with this arrival point. A special program - the Atmospheric Entry Program 
or (AEP) - was used to model the trajectory dynamics during descent and to target to the 
desired landing site. 8 The program SEPV of the DPTRAJ software set was used to target 
the interplanetary trajectory to the required atmospheric interface geometry. 9 The B - 
plane coordinates and arrival time associated with the desired descent and landing 
conditions were obtained through iterative application of these two programs. 10 SEPV 
was also used to calculate the required velocity change (AV) vector for each maneuver. 
The current best estimate of the orbit, the nominal time of the TCM and the desired TCM 
target were used together to calculate the required instantaneous AV. The times and 
overall design AV vectors for each TCM are listed in Table 4.* 


Table 3 B-PLANE AND TIME OF CLOSEST APPROACH TARGETS AND SHIFTS 

FOR TCMS 1-4 


Maneuver 

Date 

Target B*R 
(km) 

Target B*T 
(km) 

Target TCA 
(UTC on 7/4/97) 

TCM 1 

Launch + 37 days 
1/1 0/97 

-1856.93 

-4263.91 

16:53:43.20 

TCM 2 

Launch + 61 days 
2/3/97 

-2004.00 

-4554.00 

16:53:43.26 

TCM 3 

Mars - 58 days 
5/7/97 

-1855.60 

-4264.16 

16:54:22.128 

TCM 4 

Mars - 9 days 
6/25/97 

-1855.60 

-4264.16 

16:54:22.128 


Maneuver Execution Modes 


The Pathfinder propulsion system had two sets of thrusters: the axial set producing 
thrusts aligned with the spin axis of the spacecraft and the lateral set producing thrusts 
roughly perpendicular to the spin axis. Depending on the spacecraft’s attitude relative to 
the direction of the desired AV vector, one or both of these sets could be fired to add 
vectorially to the total AV. 2,7 TCM 1 was implemented using only a single session of axial 
thruster firing. This could also have been done for TCM 2, but both lateral and axial 
thrusting would be necessary for TCMs 3 and 4 given the preferred spacecraft attitude 
for late cruise. Since these TCMs were critical for Mars atmospheric entry targeting, i t 
was decided to test lateral thrusting at TCM 2. The spacecraft orientation for TCM 2 was 
chosen such that the lateral segment would match the expected magnitude of TCM 3. The 
axial segment was performed first followed by the lateral pulsing segment. The lateral 
portion of TCM 2 executed as desired, validating its use in combination with axial 
thrusting for TCMs 3 and 4. 


The maneuver design calculations assumed instantaneous application of the total AV, but the spacecraft’s 
propulsion system actually took a finite amount of time to implement these burns. For this reason, some of 
the larger maneuvers were modeled as burns of finite, non-zero duration in the OD filter, as previously 
described. 




Table 4 Design & Estimated Values of Overall Velocity Change (AV) 

Vectors for TCMs 1-4 


Maneuver 

Design 

Execution Time 
(UTC) 

Design AV Vector in EME2000 (Inertial) Coordinates 

X 

(m/sec) 

Y 

(m/sec) 

Z 

(m/sec) 

Magnitude 

(m/sec) 

TCM 1 

1/10/97 

02:00 

17.301 

-25.321 

5.899 

31.229 

TCM 2 

2/3/97 

23:00 

1.1193 

-1.1231 

0.0352 

1 .5860 

TCM 3 

5/7/97 

01 :00 

0.0762 

0.0630 

0.0362 

0.1053 

TCM 4 

6/25/97 

17:00 

0.0105 

-0.0138 

-0.0068 

0.0186 

Maneuver 

Design 

Estimated AV Vector in 

EME2000 (Inertial) Coordinates 


Execution Time 

X 

Y 

Z 

Magnitude 


(UTC) 

(m/sec) 

(m/sec) 

(m/sec) 

(m/sec) 

TCM 1 

1/10/97 

02:00 

16.6367 

-24.41 14 

5.6403 

30.077 

TCM 2 

2/3/97 

23:00 

1.1223 

-1 .1340 

0.0365 

1 .5959 

TCM 3 

5/7/97 

01 :00 

0.0764 

0.0725 

0.0382 

0.1 120 

TCM 4 

6/25/97 

17:00 

0.0105 

-0.1374 

-0.0068 

0.0186 


Contingency plans for a possible fifth TCM to be executed in the last day before entry 
were developed between TCM 2 and 3. The potential trajectory dispersions to be 
corrected by this maneuver required rather large lateral A V magnitudes. Again, it was 
desired to test performance of this lateral thrusting before any actual TCM 5 execution. 
So it was decided to perform TCM 3 in 3 segments - two lateral and one axial. A lateral 
segment was performed first to test the TCM 5 scenario. Then the axial and second lateral 
segments were performed to achieve the required total AV for TCM 3. No problems were 
encountered with the first lateral segment, verifying the design approach for TCM 5. 

TCM 4 was executed with a single lateral segment followed by an axial segment. 
Although the required AV turned out to be very small, it was decided to perform the 
maneuver in order to minimize the probability of having to do TCM 5. (In fact, TCM 5 
was not performed as explained in the section on EDL activities.) 

Maneuver Reconstruction 

The actual AVs provided by the spacecraft’s propulsion system were estimated as 
part of the on-going orbit determination process using tracking data obtained during and 
after each TCM. Table 4 shows the best estimates of the actual AV vectors for TCMs 1 - 4 . 
A comparison of these vectors with the design values reveals that maneuver execution 
was quite accurate for Mars Pathfinder. The magnitude errors achieved in flight ranged 
from 0.2% to 4%. Pointing errors, computed as the angle between the design direction 
and the best estimated direction for the TCMs, were typically less than 0.15 deg with a 
slightly larger value of just over 0.5 deg for the second lateral segment of TCM 3. Both 
of these types of errors were well within the performance requirements levied on the 
spacecraft’s propulsion system. 

The total actual AV used for Mars Pathfinder TCMs was 32.733 m/sec. This is 
slightly less than the total design value of 33.944 m/sec, primarily due to TCM 1 
executing at slightly lower than the design magnitude. The actual TCM AV expenditure 
was much smaller than the total capacity of the propulsion system and much smaller 
than the pre-launch TCM AV budget of 130-145 m/sec (99%). 2,7 This is due to the 



































































assumption of constraints for TCM 1 execution that would, in some cases, cause the 
maneuver to be executed as two separate, large axial burns. This execution mode was 
inefficient and elevated the pre-launch AV statistics. The actual TCM 1 AV direction did 
not violate these constraints, allowing the maneuver to be executed as one turn followed 
by a single axial burn. 

NAVIGATION FOR ENTRY, DESCENT AND LANDING (EDL) 

Pathfinder’s direct entry into the Martian atmosphere from its interplanetary 
trajectory led to a unique coupling between navigation knowledge, targeting to the 
atmospheric entry interface, and spacecraft activities leading up to and during the 
descent itself. The development of operational plans and procedures requiring navigation 
participation for this critical time period is discussed. A history of actual orbit 
determination results and their interaction with flight team activities at key decision 
times is presented. The final, best estimate for the spacecraft’s entry state based on all 
tracking data obtained prior to landing is given and compared to solutions available at the 
earlier decision points. This is also compared to post-landing solutions for the lander’s 
location generated from imaging data and tracking data taken after surface deployment. 

Navigation Activities and Timeline for EDL 

Pre-launch analyses had shown that no significant improvement in OD knowledge 
would occur after TCM 4 execution until the last 48 hours before entry when Mars’ 
gravitational signature would begin to be sensed in the Doppler data at levels greater 
than typical Doppler noise. Trajectory uncertainties would decrease steadily from this 
time, with the largest decreases coming in the last 24 hours. This improved trajectory 
knowledge was put to several important uses. First, it was used to update parameters i n 
the flight software controlling the mechanisms and events for a safe descent and initial 
surface deployment. These included parameters used in the primary and backup 
parachute deployment algorithms and the latitude and longitude of the surface impact 
location. Second, the predicts at the DSN station were updated with the new trajectory 
profile to increase the chances of remaining in contact with the spacecraft for as long as 
possible during descent. Finally, two opportunities were inserted to correct any 
unexpected trajectory deviations by performing a fifth maneuver -TCM 5. 

All of these considerations were factored into the operational timeline for the last 2 
days before entry shown in Figure 6. Four times were chosen as opportunities for the 
EDL parameter updates, at roughly 37 hours, 22 hours, 9 hours, and 4 hours before 
entry. Two windows were identified for possible TCM 5 execution, one at roughly 1 0 
hours (TCM 5A) and the other at roughly 6 hours before entry (TCM 5B). A single time 
was chosen for the DSN predicts update at 2.5 hours before entry (not shown in Figure 
6). The navigation team performed new orbit determination solutions at each of these 
times to update the trajectory and landing site knowledge. In the case of the EDL 
parameter updates, the team also calculated the new flight software parameter values. 
For TCM 5, the team assessed the deviations from the desired landing site and made a 
recommendation on maneuver execution. 



July 2D I July 3D I July 4D 

6 p.m. 12 a.m. 6 a.m. 12 p.m. 6 p.m. 12 a.m. 6 a.m. 12 p.m. 



□ Atmospheric flight analysis^ 

♦ Command approval meeting 

Figure 6 Timeline of Navigation Activities for EDL 

The normal maneuver design process for Pathfinder required 3-5 days of 
preparation before actually performing the burn. Obviously, there was not enough time 
to do such a customized maneuver design for the TCM 5 opportunities. Instead, a series of 
fixed AVs were generated to correct for a specific range of deviations from the desired 
landing site. It was desired to land within a surface ellipse, or footprint, 200 x 100 km 
in size with the longer axis oriented along the direction of travel of the spacecraft, or 
downtrack direction. Offsets from the nominal B-plane aimpoint along the B-vector 
mapped into downtrack dispersions from the target landing site on the surface, while 
offsets normal to the B-vector mapped into dispersions normal to the spacecraft’s 
direction of travel, or crosstrack dispersions. The B-plane ellipse representing a 6 a 
delivery accuracy for TCM 4 mapped into downtrack offsets well outside the desired 
footprint, while crosstrack dispersions remained inside it. Therefore, TCM 5 AVs were 
initially designed to cover 4 regions for both positive and negative downtrack offsets, 
giving 8 possible AV corrections for each of the two TCM 5 windows. Two corrections for 
crosstrackdeviations, or deviations normal to the B-vector, were later added for the 
TCM 5A window, bringing the total AV choices to 10 for TCM 5A and 8 for TCM 5B. 
Selection of the required AV was defined by the downtrack and crosstrack offsets or, 
equivalently the offsets along and normal to the B-vector, predicted from the latest CD 
solution. 3 

The set of overall AV s for TCM 5 downtrack and crosstrack corrections would be 
implemented using only lateral thrusting given the spacecraft’s entry attitude. Attitude 
control concerns limited the total magnitude of lateral thrusting that could be 
implemented in a single segment. TCM 5 was to be implemented as a series of one or 
more lateral AV pulsing sessions, each applying a uniform AV increment in the required 
direction. 2 to 5 of these sessions would be needed to build up the total AV for each of the 
regions of downtrack or crosstrack offsets. The sequences of commands required to 
implement each of the AV choices were designed, tested, and loaded on the spacecraft well 
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before the July 4th arrival date. This reduced the necessary activities to simply 
generating a new OD solution and identifying which of the available AVs, if any, to apply. 

The strategy for generating trajectory files for updates to the DSN predict products 
was similar to that of TCM 5 design. A series of descent trajectory files and the 
corresponding predict product sets covering the expected range of flight-path angle 
dispersions were generated several weeks in advance of landing day. A total of 21 cases 
were included, representing flight-path angles from -15.2° to -13.2° in increments of 
0.1°. The predicted entry flight-path angle from the OD solution at the designated time 
was used to select the predict set corresponding to the nearest flight-path angle value. 
Telecom equipment at the DSN stations was then reconfigured to use this set. The update 
time was chosen to allow sufficient time for this reconfiguration before receipt of the 
first signals after entry. 11 

Updated values for parameters used in the primary and backup parachute deployment 
algorithms in the EDL flight software were determined from the deceleration history 
during the initial descent into the atmosphere. The deceleration information was derived 
from propagation of sample entry states representing 3a dispersions from the predicted 
entry state of the latest OD solution. The resulting sample descent trajectories were 
generated using the NASA-Langley 6DOF POST program, 12 with the JPL 3DOF AEP 
program 8 serving as a backup and validation for the POST results. Updated values for the 
landing site latitude and longitude were taken from the AEP and POST propagations of the 
predicted entry state through the complete descent ending at surface impact. As was the 
case for TCM 5, the entry analyses associated with generating these parameter updates 
were normally quite labor-intensive. Processes for running POST and AEP to calculate 
the deceleration history and associated quantities were largely automated prior to 
landing day to provide the quick turn-around needed to meet the update schedule. Criteria 
were established for determining when the new parameter values were sufficiently 
different from current values loaded in the flight software to justify performing the 
update. These criteria along with further details on the parachute deployment algorithms 
and the entry analysis tools are given in (Ref 13). 

Orbit Determination Results for the Last 24 Hours 

The orbit determination process was the central focus for the navigation team in the 
days following execution of TCM 4 and leading up to EDL. Key decisions rested on the 
latest predictions of spacecraft entry state and landing site along with their associated 
uncertainties. Numerical results from the actual OD runs performed from July 2 to 4 , 
1997 are presented in Tables 5 and 6. Table 5 gives the predicted B-plane positions and 
times of closest approach and their associated uncertainties, as was done in Table 2. f 
Because of the importance of entry flight-path angle (FPA) at this stage of the mission, 
predictions and uncertainties for it are also given in Table 5. For a landing mission like 
Pathfinder, navigation uncertainties are more directly viewed in terms of surface 
ellipses, or footprints, centered around a predicted landing site. Centers and sizes for the 
footprints corresponding to the B-plane error ellipses are given in Table 6. A graphical 


f The times of closest approach are a convenient representation of a future interplanetary trajectory that 
assumes only gravitational forces from Mars. The spacecraft entered the Mars atmosphere, causing the 
actual trajectory to deviate from this predicted trajectory, before this closest approach time was reached. 



history of these runs is shown in the B-plane in Figure 7 and on the surface in Figure 8. 
These figures clearly illustrate the consistency of the OD results obtained at each of the 
EDL decision points. 

The largest ellipse and footprint in Figures 7 and 8 are from the run generated for 
the first EDL parameter update on July 2. This was the first time that the navigation 
team had generated results based on the final solar pressure model. Most of the new 
parameter values computed from this entry state were acceptably close to the nominal 
values, but the backup parachute deployment time differed from the nominal value by 
almost 9 sec. This was enough to justify sending the new parameter set to the spacecraft. 
Due to the consistency of the subsequent OD solutions, new parameter values remained 
very close to the values for Update #1 and no more updates were needed. The parameter 
values used for EDL Update #1 became the values used onboard the spacecraft during its 
descent to the surface; these values are listed in Table 2 of (Ref. 13). 

By the times of the TCM 5A and 5B decisions, the OD solution footprints had become 
almost completely contained in the target footprint. The default criteria indicated that 
TCM 5 was not necessary and the NAV team recommended that no TCM be performed. TCM 
4 became the final maneuver performed by the Pathfinder spacecraft. Predictions of the 
entry flight-path angle continued to yield values around -13.9° for OD solutions leading 
up to the selection time for the DSN predicts update. The trajectory file for a -13.9° 
entry FPA was chosen as the closest match given the 0.1° increments separating these 
pre-generated files. 

The NAV team continued to generate OD solutions as new tracking data were received, 
right up to the last minutes before entry. The final best estimate from pre-entry 
tracking data used 2-way Doppler data taken up to 1.3 hours before entry and range 
collected until 3.37 hours before entry. This case is labeled “All data” in Tables 5 and 6. 
The size of the landing footprint had diminished considerably with incorporation of the 
later data as shown in Figure 9. The new landing footprint was only 15x8 km assuming 
no dispersions due to atmosphere or vehicle aerodynamic mismodeling. This assumption 
was valid for earlier OD solutions where entry state uncertainties dominated the 
footprint size. As knowledge of the entry state improved, the relative significance of 
these dispersions increased. The second ellipse shown in Figure 9 represents the results 
of a Monte-Carlo analysis of the entry trajectory that incorporated appropriate 
uncertainties in aerodynamic and atmospheric conditions. The size of the predicted 
landing footprint doubles when these effects are included, as seen from the entry labeled 
“All data + descent dispersions” in Table 6. 

Independent estimates of Pathfinder’s landing location have been computed using data 
taken after surface deployment. The first of these is derived by triangulation from 
surface features identified in images taken by the lander’s camera. 13 The second estimate 
is derived by processing Doppler and range measurements from Pathfinder obtained 
between July and September 1997 along with similar data from the Viking landers. 15 
These two estimates of landing site location are shown in Table 6 under the labels “Image 
Triangulation” and “Surface Tracking Data” and marked on Figure 9. These solutions are 
reasonably consistent with the final navigation results, especially when descent 
dispersions are included in the estimate. All solutions give points roughly southwest of 
and less than 30 km from the target landing site. 



Table 5 B-PLANE STATE, TIME OF CLOSEST APPROACH, ENTRY FLIGHT-PATH 
ANGLE (FPA) AND THEIR la UNCERTAINTIES FOR EDL OD SOLUTIONS 


Solution 

Data Cutoff* 
(UTC) 

B*R 

(km) 

B*T 

(km) 

TCA 

(UTC 

on 7/4/97) 

FPA 

(deg) 

TCM 4 

6/23/97 

-1852.63 

-4250.28 

16:54:22.6 

-14.8408 

Update #1 

7/2/97 

19:46 

-1842.67 

-4275.09 

16:54:22.6 

-13.9423 

Update #2 

7/3/97 

18:48 

-1841 .28 

-4276.58 

16:54:22.6 

-13.9019 

TCM 5A 

7/4/97 

02:54 

-1844.19 

-4275.45 

16:54:22.6 

-13.8960 

Update #3 

7/4/97 

04:46 

-1844.27 

-4275.43 

16:54:22.594 

-13.8955 

TCM 5B 

7/4/97 

07:46 

-1844.52 

-4275.27 

16:54:22.589 

-13.8979 

Update #4 

7/4/97 

10:46 

-1844.58 

-4274.46 

16:54:22.516 

-13.9139 

All data 

7/4/97 

15:41 

-1850.17 

-4269.23 

16:54:22.386 

-14.0614 

Solution 

la 

SMAA 

(km) 

la 

SMIA 

(km) 

e 

(deg) 

la 

TCA 

(sec) 

la 

FPA 

(deg) 

TCM 4 

17.44 

4.192 

109.59 

0.38 

0.212 

Update #1 

13.303 

2.615 

108.016 

0.227 

0.146 

Update #2 

8.542 

2.127 

1 14.077 

0.158 

0.104 

TCM 5A 

5.816 

1 .357 

126.059 

0.1 15 

0.0899 

Update #3 

5.674 

1 .254 

127.266 

0.1 13 

0.0893 

TCM 5B 

5.380 

1.117 

129.268 

0.108 

0.0891 

Update #4 

5.106 

1 .052 

130.576 

0.103 

0.0886 

All data 

3.246 

0.176 

1 14.972 

0.023 

0.00960 


* These are the times that the last data were received on Earth; the measurement was 
actually collected at the spacecraft 10 minutes and 40 seconds earlier, which was the 
one-way light time value at entry on July 4. 


CONCLUSIONS 

All of the requirements placed on Mars Pathfinder navigation were met with ample 
margins. The lander came to rest within 30 km of the target location, well within the 
required footprint of 300 x 100 km and the smaller goal of 200 x 100 km. The final 
estimated flight-path angle at entry was -14.06°, just 0.14° off from the nominal value 
of -14.2° and well within the required 1° dispersion from the nominal. The total AV 
used for trajectory corrections was 32.7 m/sec, well under the 105 m/sec limit. 

















































































































Table 6 LANDING SITE AND 3a SURFACE FOOTPRINTS FOR EDL OD 

SOLUTIONS 


Solution 

E 

Longitude 

(deg) 

Aerocentric 
N Latitude 

(deg) 

3a 

Downtrack 

(km) 

3a 

Crosstrack 

(km) 

Footprint 

Orientation 

Angle 

(deg) 

Target 

326.9 

19.244 

N/A 

N/A 

N/A 

TCM 4 

328.5647 

19.8456 

185.661 

72.927 

10.683 

Update #1 

326.2294 

18.8238 

154.209 

51 .784 

8.21 

Update #2 

326.1944 

18.7629 

1 12.565 

37.400 

20.586 

TCM 5A 

326.1634 

18.7924 

98.301 

18.606 

29.460 

Update #3 

326.1617 

18.7929 

97.759 

16.931 

29.896 

TCM 5B 

326.1671 

18.7985 

97.600 

12.988 

30.423 

Update #4 

326.2063 

18.8291 

96.607 

12.988 

29.836 

All data 

326.5875 

19.0438 

15.116 

8.395 

84.457 

All data + descent 
dispersions 

326.48 

19.004 

41 .458 

14.807 

18.631 

Image Triangulation 

326.45 

19.14 

N/A 

N/A 

N/A 

Surface Tracking 

Data 

326.4762 

19.095 

N/A 

N/A 

N/A 
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Figure 9 Best Estimates of Pathfinder’s Landing Site 

One of the major accomplishments of Mars Pathfinder navigation was the successful 
implementation of the “enhanced filter” configuration and the precision ranging 
technique for spacecraft orbit determination. These methods allowed the necessary CD 
accuracy to be achieved using only Doppler and ranging without the addition of other data 
types such as optical images or AVLBI. Pathfinder has validated the use of this filter 
configuration for actual flight operations. 

Another significant accomplishment was the integration of existing navigation 
software for modeling an interplanetary trajectory with new software for modeling the 
descent trajectory once in the atmosphere. The JPL AEP and the NASA/Langley POST 
programs were critical additions to the navigation software suite. Coordinated used of 
these programs was necessary to target the late cruise maneuvers for a landing at the 
desired surface location. This enhanced software tool set also allowed optimization of the 
spacecraft’s descent conditions based on rapid and frequent orbit determination updates 
using the latest available tracking data in the final hours before entry. 
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